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I. INTRODUCTION 

The hallmark of the BCS model of superconductiv- 
ity in metals is the presence of an excitation gap A. 
This gap is caused by the formation of Cooper pairs de- 
scribing correlated electron pairs in time-reversed states. 
Thus, pairing correlations in superconductors tend to 
minimize the total spin of the electron system. Ferro- 
magnetic correlations, on the other hand, prefer to max- 
imize the total spin and form a macroscopic magnetic 
moment. Early worki^i 4 * 5 - predicted a state in which 
both pairing and ferromagnetic order are present, if fer- 
romagnetism is caused by localized paramagnetic impu- 
rities. The experimental observation that both states 
of matter can coexist in heavy fermion systems 6 -! 7 -! 8 - and 
high-Tc superconductors*^ came as a surprise and led 
to the search for new theoretical models to describe this 
coexistence. A BCS-like model of s-wave pairing com- 
bined with a simple Stoner-like model of ferromagnetism 
was used to derive such an intermediate state within a 
mean-field approximation?^ However, it was argued that 
such a state is unstable in the bulki 12 ' 13 ' 14 Furthermore, 
it was shown that a proper Hartree-Fock mean-field the- 
ory of the model does not support coexistence of s-wave 
superconductivity and ferromagnetism . 15 ' 16 

A similar model of BCS-like pairing and exchange in- 
teraction was shown to be valid in small metallic grains 
in the mesoscopic regime for a Thouless energy Et that is 
large compared with the single-particle mean-level spac- 
ing (5 i 17 ' 18 i 19 In such a finite-size system, a partially 
paired state with finite spin polarization exists within a 
narrow parameter regime. 20 Since this coexistence regime 
is relatively small, it would be difficult to observe it ex- 
perimentally. It has been suggested that the probability 
of spin polarization in the presence of pairing correla- 
tion may be enhanced by mesoscopic fluctuations 21 - or by 
an asymmetric spin-dependent bandwidth of the single- 
particle spectrum^ - 



Here we study the competition between ferromagnetic 
and pairing correlations in metallic grains in the crossover 
regime from a few-electron system (A <C 8) to the bulk 
(A 3> S). We use Richardson's solution of the BCS- 
like interaction 2 ^ and the known solution of the exchange 
model — to determine the ground state of the grain. For 
sufficiently small grains, there is a regime in pairing gap 
A/S and exchange coupling J s /S, in which the ground 
state is partly paired and partly polarized. We show that, 
in the presence of a Zeeman field, the exchange coupling 
at which the crossover from a pure superconducting state 
to the coexistence regime takes place decreases to values 
that can be realized in several metals. The onset of mag- 
netization with increasing exchange coupling at a given 
pairing gap corresponds to a spin jump AS > 1, followed 
by successive spin increments of AS" = 1. The magni- 
tude of the initial spin jump depends on the value of A/S. 
Similar spin jumps were found in the crossover from a su- 
perconducting state to a paramagnetic state ; 24 ' 25 where 
they are reminiscent of a first-order transition in the bulk. 
We apply a mean-field theory similar to the one used in 
Ref. [25| and compare with the exact results. In contrast 
to the exact solution, we find that the mean-field ap- 
proximation cannot describe pairing correlations in the 
intermediate regime of partial spin polarization. 

II. MODEL 

An isolated metallic grain in which the single-particle 
dynamics are chaotic and whose dimensionless Thouless 
conductance gr = Et/S is large (gr 3> 1), can be de- 
scribed by an effective universal Hamiltonia n 17 ' 18 ' 19 

H = e kcl a c ka - GP^P - J S S 2 + g[ i B HS z . (1) 

ka 

Here c\ a is the creation operator for an electron in the 
single-particle level efe with either spin up (a = +) or 
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spin down (a — — ). The one-body term in ([I]) describes 
the kinetic energy plus confining single-particle potential. 
The second term on the r.h.s. of Eq. ([1]) is a pairing in- 
teraction with strength G and where P^ = c\,c\_ is 
the pair creation operator. The third term in flj is an 
exchange interaction expressed in terms of the total spin 

operator S = Ylktra' c kcr Tcrcr ' Ckcr ' ( Ti are P aun matrices). 
The parameter J s is the exchange coupling constant (es- 
timated values of J s for a variety of materials were tab- 
ulated in Ref . l26l) . The inclusion of such an exchange in- 
teraction in quantum dots 2 -^ explained quantitatively the 
measured peak height and peak spacing statisticsj 23 i 28 
The last term on the r.h.s. of Eq. fl} describes the cou- 
pling of an external Zeeman field H (applied in the z 
direction) to the spin of the dot. Here g is the g- factor of 
the electrons in the grain (taken to be positive) and pb 
is the Bohr magneton. Orbital diamagnetism can be ne- 
glected for small grains^ The charging energy e 2 N 2 /2C 
(C is the capacitance of the grain) is a constant for a 
grain with a fixed number of electrons N and was omit- 
ted in the Hamiltonian ([TJ) . 

In this work, we do not consider mesoscopic effects that 
originate in the random matrix description of the single- 
particle Hamiltonian. To construct a typical phase dia- 
gram of a single grain, we consider a generic equidistant 
spectrum = kS with — N a < k < N a at half filling. 
Thus we have N = 2N Q for an even number of electrons 
(p = 0) and N = 2N Q + 1 for an odd number (p = 1). 



III. EXACT SOLUTION 

In the absence of a pairing interaction (G = 0), the 
Hamiltonian |T]) can be solved in closed formi 2 ^ The or- 
bital occupations hk — + Uk- commute with S 2 and 
are good quantum numbers. The empty («& = 0) and 
doubly occupied (n^ = 2) orbitals do not contribute to 
the total spin, so the total spin of the grain is obtained 
by coupling the singly occupied levels with spin 1/2 to 
total spin S and spin projection M. For a specific set 
B of b singly occupied levels, the total spin ranges from 
S=p/2toS = b/2 with each spin value having a de- 
generacy of db(S) — (5^/2) — (5+1+6/2)- A complete set 
of eigenstates is then given by |{«fc},7, S, M) where 7 
arc quantum number distinguishing between eigenstates 
with the same spim 23 i 29 

The pairing interaction can only scatter time-reversed 
pairs from doubly occupied to empty orbitals but does 
not affect the singly occupied levels (referred to as 
"blocked" levels). It is therefore sufficient to diagonalize 
the reduced BCS Hamiltonian J2ka e kC ka Ck a — GP^P us- 
ing the single-particle subspace U of empty and doubly 
occupied levels. This problem was solved by Richard- 
son^ 2 - The eigenenergies are given by 

in 

E m = ^ Efj, , (2) 



where are parameters that characterize the eigenstate 
and m = (N — b)/2 is the number of pairs. Richardson's 
parameters E a are found by solving the set of m coupled 
non-linear equations 

h^t^-Y.^ (M = l,...,m).(3) 

To each set of m doubly occupied levels at G — 0, there is 
a unique solution for Richardson's parameters at G ^ 0. 
We note that in the general case Richardson's equations 
depend on the seniority quantum numbers of the levels 
(the seniority is the number of electrons not coupled to 
spin zero). In our case, the levels are doubly degenerate 
and the seniority of a doubly occupied level is zero. 

The eigenstates constructed from the subset U of 
empty and doubly occupied levels have spin zero, so 
the total spin of the grain is determined by the spin- 
1/2 coupling of the singly occupied levels in B. The 
eigenstates of the full Hamiltonian (JXJ) are then given by 
\B, {E^}, 7, S, M) with energies of 

E = E m + *k - J S S(S + 1) + gp B HM . (4) 
keB 

In this work, we focus on the ground state of the grain 
as a function of the interaction couplings G and J s . To 
that end, we find the lowest energy E(S) in (|U) for a 
given spin S and then minimize with respect to S. The 
energy E(S) is found by choosing a set B of b = 2S 
singly occupied levels that are placed closest to the Fermi 
energy. We then populate these b levels with spin-down 
electrons, resulting in a good-spin state with spin 5* and 
spin projection M = — S. For a given set £>, we solved 
Richardson's equations using the method of Ref. I30I . 

The physical parameter dese ribing the pairing Hamil- 
tonian is A/ 5, where A is the bulk pairing gap and S 
the single-particle mean-level spacing. The low-energy 
spectrum of the grain (for J s = H = 0) is determined by 
the value of this parameter. We can truncate the total 
number of levels from N a to N r < N a , and renormalize G 
such that the low-energy spectrum of the grain remains 
approximately the same. For a picketfence spectrum, the 
renormalized coupling constant is given by 



* arcsinh(^)- 

Strictly speaking, this holds in the absence of an ex- 
change interaction. However, since the exchange interac- 
tion affects only the blocked levels, we expect the renor- 
malization ([5]) to hold as long as the number of blocked 
levels is small compared with the total number of levels in 
the band. The quality of this approximation depends on 
the choice for N r and was discussed in detail in Ref. [3l|. 
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IV. MEAN-FIELD APPROXIMATION 

We compare the findings based on the exact solution 
with a mean-field theory. The mean-field approach is 
based on a trial wave function of the form2£ 

- n 4- n { u f + v i s) 4+4-) i°) (e) 



1.2 



kEB 



jeu 



with the normalization condition (uj ) + (vj ) = 1. 
The wave function tps has b = 2S singly occupied levels 
with spin-down electrons (set B) chosen to be closest to 
the Fermi energy, and is of the BCS form within the 
remaining set of levels IA. The lowest state with spin S 
is found by minimizing the expectation value {tps I {H — 

(iN)\ips) with respect to the variational parameters v)j . 
Here \i is a chemical potential ensuring that the average 
number of particles is N. 

The mean-field energy at fixed spin S is given by 



E m{ (S) 



A 2 
G 



2 Y^' i { '', 

km 

Yjk - J S S (S+l)- gusHS , (7) 
keB 



where 



(k - M 



e k - fi) 2 + Aj 



(8) 



and As is a spin-dependent pairing gap. The gap param- 
eter and chemical potential are determined by solving the 
gap equation together with the particle number equation 



2 
G 

N 



E 

k&A 



1 



(S) 
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(9a) 



(9b) 



For an equidistant spectrum, the chemical potential can 
be determined by symmetry considerations and is given 
by fi = -(1 - p)S/2 for N 5 > A s . 

Here, we used the same approximations as in Ref. l25l 
and neglected a term in the energy E m f(S) which is pro- 
portional to {v^) 4 . The result J7]) is in agreement with 
the leading term of an expansion in the inverse number of 
electrons 1/NM. Comparing and (gj) with M = -S, 
we see that the exchange and Zeeman terms are treated 
exactly in this mean-field approximation. The ground- 
state spin in the mean-field approximation is found by 
minimizing E m [(S) in with respect to S. 



V. GROUND-STATE PHASE DIAGRAM 

The ground-state spin of the grain is determined by 
the competition between various terms in the universal 
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FIG. 1: Ground-state phase diagram in the Js/S-A/S plane 
for an even number of electrons (N — 50). Left panel: ex- 
act results. Right panel: mean-field approximation (see text). 
The phase diagrams show a superconducting (SC) phase and 
a ferromagnetic (FM) phase. The exact phase diagram also 
exhibits an intermediate regime (SC-FM) in which the ground 
state is partially polarized but still has pairing correlations. 
The intermediate regime in the mean-field phase diagram de- 
scribes a state that is partially polarized state but does not 
include pairing correlations. In particular, the dashed line 
separates an S = SC phase from an S — phase with no 
pairing correlations. The numbers shown in the intermediate 
regime are the spin values in the corresponding sectors. 

Hamiltonian. The one-body part (kinetic plus confining 
one-body potential) and pairing interaction favor mini- 
mal spin S — p/2 while exchange interaction and Zeeman 
field favor a maximally polarized state. We have studied 
the ground-state spin as a function of the three param- 
eters A/5,J s /5 and g^sH/S. Using the exact solution, 
we find three different phases: a superconducting phase 
where the number of pairs is maximal and S — p/2, a 
ferromagnetic phase where the system is fully polarized 
S = N/2 (all electrons are with spin down), and an in- 
termediate regime. The intermediate regime describes 
a partially polarized state S < N/2, in which b = 2S 
electrons reside in singly occupied levels closest to the 
Fermi energy and the remaining electrons are paired to 
give spin zero. 

The phase diagram in the A/S-J s /S plane of a grain 
with even number of electrons and in the absence of Zee- 
man field (H = 0) is presented in Fig. 1(a). For weak 
pairing, the superconducting and ferromagnetic phases 
are separated by an intermediate regime. The boundaries 
of this intermediate regime are described by two critical 
values jj 1 ^ and of the exchange interaction that are 
function of A/6. The critical value Js /S is a monotoni- 
cally increasing function of A/5, i.e., a stronger exchange 
interaction is required to polarize a grain with stronger 

(2) 

pairing correlations. However, J s /S is almost insensi- 
tive to A/ 6. The intermediate regime shrinks at larger 
A/ 6 and eventually disappears above A/ 5 ~ 3. For 
stronger pairing correlations, the superconducting phase 
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makes a direct transition to the ferromagnetic phase. In 
this regime (not shown in Fig. 1), the phase boundary 
exhibits a strong dependence on the bandwidth N D . 

For comparison, we show the mean-field results in 
Fig. 1(b). We observe that the mean-field results are 
qualitatively different from the exact solution. The re- 
gion to the right of the thick solid line and dashed line 
describes an superconducting phase with Ao 7^ 0. How- 
ever, there is no superconducting solution (i.e., Ao = 0) 
for A/S < 0.28. Furthermore, in each of the partially po- 
larized regions with spin < S < N/2, the correspond- 
ing pairing gap vanishes A5 = and there are no pairing 
correlations present. While solutions with A5 ^ exist, 
they occur for values of A/S for which a higher spin state 
with no pairing correlations has lower energy (because of 
the exchange interaction). For example, a solution with 
Ai ^ exists only for A/S > 2.1^ However, at this 
strength of pairing correlation we observe a direct tran- 
sition from S = to S = 4 with As — as the energy 
of the lowest S = 4 state with no pairing correlations 
is lower than the paired S = 1 state. As a result, the 
boundaries which separate different spin phases are flat, 
e.g., independent of the pairing strength. In fact, in the 
mean-field approximation the ground-state wave function 
is a Slater determinant through the whole intermediate 
regime of partial spin polarization. Thus no coexistence 
of pairing and spin polarization is observed within the 
mean-field approach. 

In contrast, the exact solution shows that pairing cor- 
relations are present as long as the system is not fully 
polarized. This can be seen in the shift of the spin tran- 
sition lines to higher values of the exchange interaction 
strength as the pairing gap A/S is increased. Thus, the 
exact solution predicts a regime in which pairing correla- 
tions and spin polarization coexist. In the following, we 
only discuss results obtained from the exact solution. 

More detailed phase diagrams for H = are shown in 
the top row of Fig. 2 for both grains with even [panel (a)] 
and odd [panel (b)] number of electrons. For weak pairing 
we observe an odd-even effect (in number of electrons). 

In particular, the critical value Js 1 ^ is larger for the odd 
grain, even though the presence of a blocked level in the 
odd superconducting phase weakens pairing correlations 
in the odd grain. This is because increasing the spin from 
1/2 to 3/2 in the odd grain costs more one-body energy 
than increasing the spin from to 1 in the even grain. 

The phase boundaries for a finite Zeeman field 
gfisH/S — 2.6 are shown in the bottom row of Fig. 2. 
The effect of a Zeeman field is twofold. First, it helps 
polarizing the grain, making the value of Js 1 ^ for a given 
pairing gap smaller. Second, at given exchange strength 
J s /S, it increases the critical value of A/S at which par- 
tial spin polarization is destroyed. Both effects together 
increase the size of the intermediate regime in the A/5— 
J s /S plane. 
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FIG. 2: Phase diagrams in the J 3 /5-A/5 plane at a fixed 
Zeeman field gg.aH = (top panels) and gg.BH/5 = 2.6 
(bottom panels) for an even grain (left panels) and for an 
odd grain (right panels). The numbers denote the spin in 
each sector. 

VI. SPIN JUMPS 

As we increase the exchange coupling constant J s /S at 
fixed A/5 and Zeeman field, the spin increases by discrete 
steps from its minimal value S — p/2 to its maximal 
value of S = N/2. In the absence of pairing (A = 0), 
the transition from spin S to spin S + 1 occurs for an 
exchange coupling of 

(2S+l)-gy B H/5 
J s /5 = 2S + 2 A = . (10) 

The corresponding A = phase diagrams in the 
gUsH/S-A/S plane are shown in Fig. 3(a) and 3(b) for 
even and odd grains, respectively. In particular, the 
phase boundaries are given by = S(p+ l)/(p+ 2) — 
gp, B H/(p + 2) and J s (2) = 5{N - l)/N - gfi B H/N. The 
ground-state spin increases as a function of J s in steps of 
AS* = 1. An interesting qualitative change in the pres- 
ence of pairing correlations is the possibility of a spin 
jump AS > 1. For A/5 < 0.6, the ground-state spin 
still increases in steps of AS* = 1 versus J s . However, for 
0.6 < A/5 < 0.8, the ground-state spin jumps from to 
2 within the range 0.87 < J s /5 < 0.9. The size of the 
first-step spin jump gets larger with increasing A/S. All 
subsequent steps are of size one [see Fig. 2(a)]. 

A similar effect was observed when superconductivity 
breaks down due to the presence of a large external Zee- 
man fields The experimental findings were qualitatively 
explained using the mean-field theory we discussed pre- 
viously (but without the inclusion of an exchange inter- 
action) . It was concluded in Ref. 25 that the first-order 
phase transition from a superconductor to a paramagnet, 
observed in thin films, is "softened" in metallic grains. 



FIG. 3: Phase diagrams in the J a / S-gusH / 8 plane at fixed 
A/5 = (top row) and A/ 5 = 2 (bottom row) for an even 
grain (left column) and for an odd grain (right column). Num- 
bers denote the spin in each sector. 



Here we have shown that spin jumps also occur in the 
presence of exchange correlations. In the absence to an 
external Zeeman field, these spin jumps are predicted to 
occur at exchange coupling values J s /S > 0.87. Such ex- 
change coupling values are significantly lar ger than the 
values for most metals (see Fig. 9 in Ref. |26I ). More- 
over, the exchange is an intrinsic material property and 
is difficult to tune experimentally. 

The regime of spin jumps can be tuned to lower and 
more typical values of J s by applying an external Zeeman 
field. We have already seen in Figs. 2(c) and 2(d) that a 
relatively weak Zeeman field increases the size of the in- 
termediate regime. It also means that spin jumps can be 
observed at smaller values of the exchange strength that 
are accessible to experiments. This is demonstrated in 
Fig. 3(c) and 3(d) where phase diagrams in the g/j.BH/5- 
A/S plane are shown for a given pairing gap of A/5 = 2. 
For example, a Zeeman field of g/isH/S « 2 is sufficient 
to lower the critical exchange strength for the 0^2 spin 
jump to J s /S w 0.55 at A/S = 2 [see Fig. 3(c)] as com- 
pared to J s — 0.89(5 at A = 0.7(5 without Zeeman field 
[see Fig. 2(a)]. 

The idea of a Zeeman-field tuning of the values of ex- 
change coupling where spin jumps occur is best illus- 
trated in Fig. 4, where spin staircase functions are shown 
versus J s /S. In the presence of pairing correlations and 
the absence of Zeeman field, the ground-state spin stair- 
case is shifted to the right and compressed as A/S in- 
creases [see Fig. 4(a)], reflecting the fact that the inter- 
mediate region shrinks [see Fig. 2(a) and 2(b)]. For an 
even grain with A/S = 0.7, a spin jump of AS* = 2 sets 
in at J s w 0.895, while for A/<5 = 0.9, a spin jump of 
AS = 3 occurs at J s « 0.92(5. For a finite Zeeman field 
of gfisH/S = 2, the spin staircase functions that exhibit 



FIG. 4: Ground-state spin versus exchange coupling J„/S for 
an even grain at a fixed Zeeman field strength 2\ibH = 
(top panel) and g^BH/S = 2 (bottom panel). Solid lines 
correspond to a grain with no pairing correlations A/S = 0. 
The dashed lines describe staircase functions with a spin jump 
of AS = 2 for A/<5 = 0.7 (top) and A/S = 2 (bottom). 
The dotted lines correspond to staircases with a spin jump of 
AS = 3 for A/S = 0.9 (top) and A/<5 = 2.3 (bottom). 

similar spin jumps are shifted to smaller values of the ex- 
change strength but larger values of the pairing gap [see 
Fig. 4(b)]. Spin jumps of AS" = 2 (AS = 3) occur at 
J s /S = 0.55 (J s /8 = 0.64) and A/S = 2 (A/<5 = 2.3). 

In the relevant experimental situation of a fixed ex- 
change interaction strength, the critical value of A/S at 
which spin jumps occur as well as the size of these jumps 
increase at larger values of gusH/S. The ratio A/<5 (for 
the given metal) can be made larger by studying a larger 
grain, hence reducing the mean level spacing S. 

As an example, niobium has an exchange interaction 
strength of J s /S w 0.4 (see, e.g., Fig. 9 in Ref. H). 
Without an external Zeeman field, the ground-state spin 
will be minimal (S = p/2) at all values of A/S [see 
Fig. 2(a) and 2(b)]. At a Zeeman field g^ B H/5 = I, 
the ground-state spin of niobium changes from to I 
at A/S = 0.66. However, at g^sH/S — 2.6, a spin 
jump of AS — 2 occurs from to 2 at A/S = 2.15 
[see Fig. 2(c)]. For these two values of A/<5, we can 
roughly estimate the corresponding critical size of the 
metallic grain given the bulk gap value A = 3.05 meV 
and Fermi momentum fcp = 11.8 nm" 1 of niobium. In a 
Fermi gas model, the mean-level spacing is related to the 
volume of the grain by S = 2ir 2 h 2 / '(rnk-pV). Assuming 
a hemispheric grain with radius r, we have the relation 

1/3 

rNb ~ 2.7nm (ANb/<5) • Thus, the Hamiltonian (fT|) 
with an equidistant spectrum predicts a — > 1 spin tran- 
sition for a niobium grain of radius r w 2.35 nm and Zee- 
man field of gfisH — 4.62 meV, and a — > 2 spin jump 
at r w 3.48 nm and a Zeeman field of g^sH = 3.69 meV. 
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VII. CONCLUSION 

We have shown that there exists a small region in the 
ground-state phase diagram of a small metallic grain in 
which pairing correlations and ferromagnetism coexist. 
This coexistence regime becomes larger (in the J s /5-A/5 
plane) and therefore more accessible to experiments in 
the presence of a finite Zeeman field. In particular, we 
propose that for a given exchange constant (determined 
by the material used), spin jumps can be observed by tun- 
ing a Zeeman field. We have also shown that a quantita- 
tive study of the intermediate regime requires the use of 
the exact solution. Furthermore, the mean-field approx- 
imation is qualitatively different in that it does not pre- 
dict any pairing correlations in the intermediate regime 



of partial spin polarization. 

In this work, we have ignored mesoscopic fluctuations 
and focused on a grain with an equidistant single-particle 
spectrum. It would be interesting to study how meso- 
scopic fluctuations affect the boundaries of the interme- 
diate phase and the size of spin jumps. 
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